Nuclear and mitochondrial genomes of Polylopha cassiicola: the first assembly in Chlidanotinae (Tortricidae).

Tortricidae is one of the largest families in Lepidoptera, including subfamilies of Tortricinae, Olethreutinae, and Chlidanotinae. Here, we assembled the gap-free genome for the subfamily Chlidanotinae using Illumina, Nanopore, and Hi-C sequencing from Polylopha cassiicola, a pest of camphor trees in southern China. The nuclear genome is 302.03 Mb in size, with 36.82% of repeats and 98.4% of BUCSO completeness. The karyotype is 2n = 44 for males. We identified 15412 protein-coding genes, 1052 tRNAs, and 67 rRNAs. We also determined the mitochondrial genome of this species and annotated 13 protein-coding genes, 22 tRNAs, and one rRNA. These high-quality genomes provide valuable information for studying phylogeny, karyotypic evolution, and adaptive evolution of tortricid moths.

A chromosome-level genome for the flower thrips Frankliniella intonsa.

The flower thrips Frankliniella intonsa (Thysanoptera: Thripidae) is a common insect found in flowers of many plants. Sometimes, F. intonsa causes damage to crops through direct feeding and transmission of plant viruses. Here, we assembled a chromosomal level genome of F. intonsa using the Illumina, Oxford Nanopore (ONT), and Hi-C technologies. The assembled genome had a size of 209.09 Mb, with a contig N50 of 997 bp, scaffold N50 of 13.415 Mb, and BUSCO completeness of 92.5%. The assembled contigs were anchored on 15 chromosomes. A set of 14,109 protein-coding genes were annotated in the genome with a BUSCO completeness of 95.0%. The genome contained 491 non-coding RNA and 0.57% of interspersed repeats. This high-quality genome provides a valuable resource for understanding the ecology, genetics, and evolution of F. intonsa, as well as for controlling thrips pests.

Plasticity of cold and heat stress tolerance induced by hardening and acclimation in the melon thrips.

Extreme temperatures threaten species under climate change and can limit range expansions. Many species cope with changing environments through plastic changes. This study tested phenotypic changes in heat and cold tolerance under hardening and acclimation in the melon thrips, Thrips palmi Karny (Thysanoptera: Thripidae), an agricultural pest of many vegetables. We first measured the critical thermal maximum (CTmax) of the species by the knockdown time under static temperatures and found support for an injury accumulation model of heat stress. The inferred knockdown time at 39 °C was 82.22 min. Rapid heat hardening for 1 h at 35 °C slightly increased CTmax by 1.04 min but decreased it following exposure to 31 °C by 3.46 min and 39 °C by 6.78 min. Heat acclimation for 2 and 4 days significantly increased CTmax at 35 °C by 1.83, and 6.83 min, respectively. Rapid cold hardening at 0 °C and 4 °C for 2 h, and cold acclimation at 10 °C for 3 days also significantly increased cold tolerance by 6.09, 5.82, and 2.00 min, respectively, while cold hardening at 8 °C for 2 h and acclimation at 4 °C and 10 °C for 5 days did not change cold stress tolerance. Mortality at 4 °C for 3 and 5 days reached 24.07 % and 43.22 % respectively. Our study showed plasticity for heat and cold stress tolerance in T. palmi, but the thermal and temporal space for heat stress induction is narrower than for cold stress induction.

The complete mitochondrial genome of Aphidius colemani (Hymenoptera: Braconidae: Aphidiinae).

The genome-level features are crucial genetic resources for species identification and phylogenetic analysis. Here, the complete mitochondrial genome of Aphidius colemani Viereck 1912 (Hymenoptera: Braconidae: Aphidiinae) was sequenced, determined and analyzed. The circular genome is 16,372 bp in length with an overall base composition of 38.9% for A, 46.2% for T, 6.7% for C, and 8.2% for G. The mitochondrial genome of A. colemani contained 13 protein-coding genes that initiated by the ATN codon, 22 transfer RNA genes, two ribosomal RNA genes (rRNAs), and a control region (CR). It shared the same gene arrangement patterns that occurred in two tRNA clusters of trnI-trnQ-trnM and trnW-trnC-trnY with Aphidius gifuensis. Phylogenetic analyses using Bayesian inference and Maximum-likelihood methods supported that the two species of Aphidiinae formed a clade and sister to other subfamilies of Braconidae.

Rapid and Repeated Climate Adaptation Involving Chromosome Inversions following Invasion of an Insect.

Following invasion, insects can become adapted to conditions experienced in their invasive range, but there are few studies on the speed of adaptation and its genomic basis. Here, we examine a small insect pest, Thrips palmi, following its contemporary range expansion across a sharp climate gradient from the subtropics to temperate areas. We first found a geographically associated population genetic structure and inferred a stepping-stone dispersal pattern in this pest from the open fields of southern China to greenhouse environments of northern regions, with limited gene flow after colonization. In common garden experiments, both the field and greenhouse groups exhibited clinal patterns in thermal tolerance as measured by critical thermal maximum (CTmax) closely linked with latitude and temperature variables. A selection experiment reinforced the evolutionary potential of CTmax with an estimated h2 of 6.8% for the trait. We identified 3 inversions in the genome that were closely associated with CTmax, accounting for 49.9%, 19.6%, and 8.6% of the variance in CTmax among populations. Other genomic variations in CTmax outside the inversion region were specific to certain populations but functionally conserved. These findings highlight rapid adaptation to CTmax in both open field and greenhouse populations and reiterate the importance of inversions behaving as large-effect alleles in climate adaptation.

Chromosome-level genome assembly of the Japanese sawyer beetle Monochamus alternatus.

The Japanese sawyer beetle Monochamus alternatus (Coleoptera: Cerambycidae) is a pest in pine forests and acts as a vector for the pine wood nematode Bursaphelenchus xylophilus, which causes the pine wilt disease. We assembled a high-quality genome of M. alternatus at the chromosomal level using Illumina, Nanopore, and Hi-C sequencing technologies. The assembled genome is 767.12 Mb, with a scaffold N50 of 82.0 Mb. All contigs were assembled into ten pseudo-chromosomes. The genome contains 63.95% repeat sequences. We identify 16, 284 protein-coding genes in the genome, of which 11,244 were functionally annotated. The high-quality genome of M. alternatus provides an invaluable resource for the biological, ecological, and genetic study of this beetle and opens new avenues for understanding the transmission of pine wood nematode by insect vectors.

Stable composition of gut microbiome in the Asian ladybeetle Coccinella septempunctata reared on natural and artificial diets.

The Asian ladybeetle, Coccinella septempunctata, is an important insect of predatory natural enemy, which has a strong control effect and application prospects for aphids, whiteflies, mealybugs, and other small-sized pests of agriculture and forestry crops. Gut microbiota composition posed impacts on development of insects. In order to clarify the effect of artificial feed feeding on the intestinal microbial species and structure of C. septempunctata, we compared the intestinal microbial flora of C. septempunctata reared on bean aphids and artificial food for 15 days. Results show that Proteobacteria was the dominant component in all groups at phylum level, Rhodobacter, Methylovigula, Burkholderia, and Bradyrhizobium were the dominant bacteria among all groups at genus level. As to the differences in bacterial community structure and diversity, there is no significant difference between Shannon index and Simpson index, the principal components analysis of the bacterial communities, and the samples were roughly distributed in different regions. After 15 days of feeding, artificial diet did not significantly reduce the microbial diversity of the gut of C. septempunctata compared to the aphid group, and there was no significant effect on the abundance of dominant flora in the gut of C. septempunctata, C. septempunctata gut has a similar core microbiota. This study clarifies the effects in intestinal microbial diversity and composition structure of the C. septempunctata with artificial diet, and provides a theoretical basis for understanding the intestinal microorganisms and optimizating the artificial diet of C. septempunctata.

Chromosome-level genome assembly of the giant ladybug Megalocaria dilatata.

The giant ladybug Megalocaria dilatata (Fabricius) is a potential biocontrol agent and a valuable model for coccinellid genomics and evolutionary biology. However, the lack of a reference genome for M. dilatata has impeded further explorations into its evolution and constrained its use in pest management. Here, we assembled and annotated a high-quality, chromosome-level genome of M. dilatata. The resulting assembly spans 772.3 Mb, with a scaffold N50 of 72.48 Mb and a GC content of 34.23%. The Hi-C data aided in anchoring the assembly onto 10 chromosomes ranging from 43.35 to 108.16 Mb. We identified 493.33 Mb of repeat sequences, accounting for 63.88% of the assembled genome. Our gene prediction identified 25,346 genes, with 81.89% annotated in public protein databases. The genome data will provide a valuable resource for studying the biology and evolution of Coccinellidae, aiding in pest control strategies and advancing research in the field.

Chromosome-level genome assembly of an agricultural pest Zeugodacus tau (Diptera: Tephritidae).

The fruit fly Zeugodacus tau (Diptera: Tephritidae) is a major pest of melons and other cucurbits in Southeast Asia. In this study, we used Illumina, Nanopore, and Hi-C sequencing technologies to assemble a reference genome of Z. tau at the chromosomal level. The assembled genome was 421.79 Mb and consisted of six chromosomes (one X-chromosome + five autosomes). The contig N50 was 4.23 Mb. We identified 20,922 protein-coding genes, of which 17,251 (82.45%) were functionally annotated. Additionally, we found 247 rRNAs, 435 tRNAs, 67 small nuclear RNAs, and 829 small RNAs in the genome. Repetitive elements accounted for 55.30 Mb (13.15%) of the genome. This high-quality genome assembly is valuable for evolutionary and genetic studies of Z. tau and its relative species.

Effects of fungicides on fitness and Buchnera endosymbiont density in Aphis gossypii.

BACKGROUND: Several agricultural fungicides are known to affect insect pests directly and these effects may be transgenerational and mediated through impacts on endosymbionts, providing opportunities for pest control. The cotton aphid Aphis gossypii is a polyphagous pest that can cause large crop yield losses. Here, we tested the effects of three fungicides, pyraclostrobin, trifloxystrobin and chlorothalonil, on the fitness and Buchnera endosymbiont of A. gossypii. RESULTS: The formulations of trifloxystrobin and pyraclostrobin, and the active ingredient of pyraclostrobin produced dose-dependent mortality in A. gossypii, whereas there was no dose-dependent mortality for chlorothalonil. The formulations of trifloxystrobin and pyraclostrobin significantly reduced the lifespan and fecundity of A. gossypii, and increased the density of Buchnera in the parental generation but not the (unexposed) F1 . When the active ingredient of pyraclostrobin was tested, the lifespan of the F0 generation was also reduced, but the density of Buchnera was not, indicating that non-insecticidal chemicals in the fungicide formulation may affect the density of the endosymbiont of A. gossypii. There was no transgenerational effect of the active ingredient of pyraclostrobin on the lifespan and Buchnera of (unexposed) F1 . CONCLUSIONS: Our results suggest that formulations of two strobilurin fungicides have immediate impacts on the fitness of A. gossypii, and chemicals in the formulation impact the density of the primary Buchnera endosymbiont. Our study highlights the potential effects of non-insecticidal chemicals of fungicides on aphid pests and their primary endosymbionts but direct connections between fitness and Buchnera densities remain unclear. © 2023 Society of Chemical Industry.

Testing for adaptive changes linked to range expansion following a single introduction of the fall webworm.

Adaptive evolution following colonization can affect the impact of invasive species. The fall webworm (FWW) invaded China 40 years ago through a single introduction event involving a severe bottleneck and subsequently diverged into two genetic groups. The well-recorded invasion history of FWW, coupled with a clear pattern of genetic divergence, provides an opportunity to investigate whether there is any sign of adaptive evolution following the invasion. Based on genome-wide SNPs, we identified genetically separated western and eastern groups of FWW and correlated spatial variation in SNPs with geographical and climatic factors. Geographical factors explained a similar proportion of the genetic variation across all populations compared with climatic factors. However, when the two population groups were analysed separately, environmental factors explained more variation than geographical factors. SNP outliers in populations of the western group had relatively stronger response to precipitation than temperature-related variables. Functional annotation of SNP outliers identified genes associated with insect cuticle protein potentially related to desiccation adaptation in the western group and genes associated with lipase biosynthesis potentially related to temperature adaptation in the eastern group. Our study suggests that invasive species may maintain the evolutionary potential to adapt to heterogeneous environments despite a single invasion event. The molecular data suggest that quantitative trait comparisons across environments would be worthwhile.

Interspecific and intraspecific variation in susceptibility of two co-occurring pest thrips, Frankliniella occidentalis and Thrips palmi, to nine insecticides.

BACKGROUND: Field control of pest thrips mainly relies on insecticides, but the toxicity of insecticides can vary among thrips species and populations. In this study, we examined the susceptibility of multiple field populations of two thrips pests, Frankliniella occidentalis, and Thrips palmi, that often co-occur on vegetables, to nine insecticides belonging to seven subgroups. RESULTS: The highest level of variation in susceptibility among F. occidentalis populations was for spinetoram (73.92 fold difference between most resistant and most susceptible population), followed by three neonicotinoids (8.06-15.99 fold), while among T. palmi populations, it was also for spinetoram (257.19 fold), followed by emamectin benzoate, sulfoxaflor, and acetamiprid (23.64-45.50 fold). These findings suggest evolved resistance to these insecticides in some populations of the two thrips. One population of F. occidentalis had a particularly high level of resistance overall, being the most resistant for five of the nine insecticides tested. Likewise, a population of T. palmi had high resistance to all nine insecticides, again suggesting the evolution of resistance to multiple chemicals. For F. occidentalis, the LC95 values of most populations were higher than the field-recommended dosage for all insecticides except chlorfenapyr and emamectin benzoate. For several T. palmi populations, the LC95 values also tended to be higher than recommended dosages, except in the case of emamectin benzoate and spinetoram. CONCLUSIONS: Our study found interspecific and intraspecific variations in the susceptibility of two thrips to nine insecticides and multiple resistance in some populations, highlighting the need for ongoing monitoring and resistance management. © 2023 Society of Chemical Industry.

The Complete Mitochondrial Genome of the Chinese White Wax Scale Insect, Ericerus pela Chavannes (Hemiptera: Coccidae), with Novel Gene Arrangement and Truncated tRNA Genes.

The Chinese white wax scale insect, Ericerus pela Chavannes (Hemiptera: Coccidae), is one of the scale insects with great economic value and has been dispersed and reared in China for over one thousand years. Its mitochondrial genome provides essential information for the molecular identification and genetic study of this species. We assembled the complete mitochondrial genome of E. pela based on PacBio sequencing and analyzed its genomic features. The genome was 17,766 bp in length with 13 protein-coding genes, 22 tRNAs, and two rRNA genes. The analysis results showed E. pela had significant gene rearrangements involving tRNAs compared with other Coccoidea species. Furthermore, E. pela's nine tRNAs were identified to have obvious truncated structures. The phylogenetic tree compiled of the species showed a long branch of the Coccoidea lineage, which indicated the high evolutionary rate in this group. Our study revealed the mitochondrial characteristics of E. pela and enriched the mitochondrial genetic information on Coccoidea species. It also determined the occurrence of gene rearrangement for the species in this superfamily.

A high-throughput KASP assay provides insights into the evolution of multiple resistant mutations in populations of the two-spotted spider mite Tetranychus urticae across China.

BACKGROUND: The two-spotted spider mite (TSSM), Tetranychus urticae (Acari: Tetranychidae), is a cosmopolitan phytophagous pest in agriculture and horticulture. It has developed resistance to many acaricides by target-site mutations. Understanding the status and evolution of resistant mutations in the field is essential for resistance management. Here, we applied a high-throughput Kompetitive allele-specific polymerase chain reaction (KASP) method for detecting six mutations conferring resistance to four acaricides of the TSSM. We genotyped 3274 female adults of TSSM from 43 populations collected across China in 2017, 2020, and 2021. RESULTS: The KASP genotyping of 24 testing individuals showed 99% agreement with Sanger sequencing results. KASP assays showed that most populations had a high frequency of mutations conferring avermectin (G314D and G326E) and pyridaben (H92R) resistance. The frequency of mutation conferring bifenazate (A269V and G126S) and etoxazole (I1017F) resistance was relatively low. Multiple mutations were common in the TSSM, with 70.2% and 24.6% of individuals having 2-6 and 7-10 of 10 possible resistant alleles, respectively. No loci were linked in most populations among the six mutations, indicating the development of multiple resistance is mainly by independent selection. However, G314D and I1017F on the nuclear genome deviated from Hardy-Weinberg equilibrium in most populations, indicating significant selective pressure on TSSM populations by acaricides or fitness cost of the mutations in the absence of acaricide selection. CONCLUSION: Our study revealed that the high frequency of TSSMs evolved multiple resistant mutations in population and individual levels by independent selection across China, alarming for managing multiple-acaricides resistance. © 2023 Society of Chemical Industry.

KASP genotyping and semi-quantitation of G275E mutation in the α6 subunit of Thrips palmi nAChR gene conferring spinetoram resistance.

BACKGROUND: Pesticide resistance is a long-standing and growing problem in the chemical control of invertebrate pests. Molecular diagnostic methods can facilitate pesticide resistance management by accurately and efficiently detecting resistant mutations and their frequency. In this study, the kompetitive allele specific PCR (KASP) approach, a technology for high-throughput single nucleotide polymorphism (SNP) genotyping, is validated as a useful method for characterizing genotypes at a pesticide-resistance locus for the first time. We focus on the spinetoram resistance mutation of G275E in the nicotinic acetylcholine receptor alpha 6 (nAChR α6) subunit gene of Thrips palmi. RESULTS: Of the 341 individuals of Thrips palmi tested, 98.24% were successfully genotyped, with 100% concordance with Sanger sequencing results. We then quantitatively mixed genomic DNA of known genotypes to establish 21 DNA mixtures with a resistant allele frequency ranging from 0 to 100% at steps of 5%. The linear discriminant analysis (LDA) showed that 75.8% of original grouped cases were correctly classified; six groups had no overlap in membership (resistant allele frequency: 0%, 5%, 10-75%, 80-85%, 90-95%, and 100%). When we chose 11 pooled samples with 10% steps for LDA, 84.4% of original grouped cases were correctly classified; seven groups had no overlap in membership (0%, 10%, 20-30%, 40-70%, 80%, 90%, 100%). The results indicated that KASP applied to pooled samples may provide a semi-quantitative estimate of resistance. CONCLUSIONS: Our study points to the suitability of KASP for high-throughput genotyping of genotypes affecting pesticide resistance and semi-quantitative assessments of resistance allele frequencies in populations. © 2023 Society of Chemical Industry.

A comprehensive assessment of insecticide resistance mutations in source and immigrant populations of the diamondback moth Plutella xylostella (L.).

BACKGROUND: The diamondback moth (DBM) Plutella xylostella has developed resistance to almost all insecticides used to control it. Populations of DBM in temperate regions mainly migrate from annual breeding areas. However, the distribution pattern of insecticide resistance of DBM within the context of long-distance migration remains unclear. RESULTS: In this study, we examined the frequency of 14 resistance mutations for 52 populations of DBM collected in 2010, 2011, 2017 and 2018 across China using a high-throughput KASP genotyping method. Mutations L1041F and T929I conferring pyrethroid resistance, and mutations G4946E and E1338D conferring chlorantraniliprole resistance were near fixation in most populations, whereas resistant alleles of F1020S, M918I, A309V and F1845Y were uncommon or absent in most populations. Resistance allele frequencies were relatively stable among different years, although the frequency of two mutations decreased. Principal component analysis based on resistant allele frequencies separated a southern population as an outlier, whereas the immigrants clustered with other populations, congruent with the migration pattern of northern immigrants coming from the Sichuan area of southwestern China. Most resistant mutations deviated from Hardy-Weinberg equilibrium due to a lower than expected frequency of heterozygotes. The deviation index of heterozygosity for resistant alleles was significantly higher than the index obtained from single nucleotide polymorphisms across the genome. These findings suggest heterogeneous selection pressures on resistant mutations. CONCLUSION: Our results provide a picture of resistant mutation patterns in DBM shaped by insecticide usage and migration of this pest, and highlight the widespread distribution of resistance alleles in DBM. © 2022 Society of Chemical Industry.

Close Kin Dyads Indicate Intergenerational Dispersal and Barriers.

AbstractThe movement of individuals through continuous space is typically constrained by dispersal ability and dispersal barriers. A range of approaches have been developed to investigate these. Kindisperse is a new approach that infers recent intergenerational dispersal (σ) from close kin dyads and appears particularly useful for investigating taxa that are difficult to observe individually. This study, focusing on the mosquito Aedes aegypti, shows how the same close kin data can also be used for barrier detection. We empirically demonstrate this new extension of the method using genome-wide sequence data from 266 Ae. aegypti. First, we use the spatial distribution of full-sib dyads collected within one generation to infer past movements of ovipositing female mosquitoes. These dyads indicated the relative barrier strengths of two roads and performed favorably against alternative genetic methods for detecting barriers. We then use Kindisperse to quantify recent intergenerational dispersal (σ=81.5-197.1 m generation-1/2) from the difference in variance between the sib and the first cousin spatial distributions and, from this, estimate effective population density (ρ=833-4,864 km-2). Dispersal estimates showed general agreement with those from mark-release-recapture studies. Barriers, σ, ρ, and neighborhood size (331-526) can inform forthcoming releases of dengue-suppressing Wolbachia bacteria into this mosquito population.

Barley Yellow Dwarf Virus Influences Its Vector's Endosymbionts but Not Its Thermotolerance.

The barley yellow dwarf virus (BYDV) of cereals is thought to substantially increase the high-temperature tolerance of its aphid vector, Rhopalosiphum padi, which may enhance its transmission efficiency. This is based on experiments with North American strains of BYDV and R. padi. Here, we independently test these by measuring the temperature tolerance, via Critical Thermal Maximum (CTmax) and knockdown time, of Australian R. padi infected with a local BYDV isolate. We further consider the interaction between BYDV transmission, the primary endosymbiont of R. padi (Buchnera aphidicola), and a transinfected secondary endosymbiont (Rickettsiella viridis) which reduces the thermotolerance of other aphid species. We failed to find an increase in tolerance to high temperatures in BYDV-infected aphids or an impact of Rickettsiella on thermotolerance. However, BYDV interacted with R. padi endosymbionts in unexpected ways, suppressing the density of Buchnera and Rickettsiella. BYDV density was also fourfold higher in Rickettsiella-infected aphids. Our findings indicate that BYDV does not necessarily increase the temperature tolerance of the aphid transmission vector to increase its transmission potential, at least for the genotype combinations tested here. The interactions between BYDV and Rickettsiella suggest new ways in which aphid endosymbionts may influence how BYDV spreads, which needs further testing in a field context.

Genome-wide and expression-profiling analyses of the cytochrome P450 genes in Tenebrionidea.

Cytochrome P450 monooxygenases (CYPs) are present in almost all areas of the tree of life. As one of the largest and most diverse superfamilies of multifunctional enzymes, they play important roles in the metabolism of xenobiotics and biosynthesis of endogenous compounds, shaping the success of insects. In this study, the CYPome (an omics term for all the CYP genes in a genome) diversification was examined in the four Tenebrionidea species through genome-wide analysis. A total of 483 CYP genes were identified, of which 103, 157, 122, and 101 were respectively deciphered from the genomes of Tebebrio molitor, Asbolus verucosus, Hycleus cichorii and Hycleus phaleratus. These CYPs were classified into four major clans (mitochondrial, CYP2, CYP3, and CYP4), and clans CYP3 and CYP4 are most diverse. Phylogenetic analysis showed that most CYPs of these Tenebrionidea beetles from each clan had a very close 1:1 orthology to each other, suggesting that they originate closely and have evolutionally conserved function. Expression analysis at different developmental stages and in various tissues showed the life stage-, gut-, salivary gland-, fat body-, Malpighian tubule-, antennae-, ovary- and testis-specific expression patterns of T. molitor CYP genes, implying their various potential roles in development, detoxification, immune response, digestion, olfaction, and reproduction. Our studies provide a platform to understand the evolution of Tenebrionidea CYP gene superfamily, and a basis for further functional investigation of the T. molitor CYPs involved in various biological processes.

Comparative genomic analysis of carboxylesterase genes in Tenebrio molitor and other four tenebrionids.

Carboxylesterases (COEs) have various functions in wide taxons of organisms. In insects, COEs are important enzymes involved in the hydrolysis of a variety of ester-containing xenobiotics, neural signal transmission, pheromone degradation, and reproductive development. Understanding the diversity of COEs is basic to illustrate their functions. In this study, we identified 53, 105, 37, and 39 COEs from the genomes of Tenebrio molitor, Asbolus verucosus, Hycleus cichorii, and H. phaleratus in the superfamily of Tenebrionidea, respectively. Phylogenetic analysis showed that 234 COEs from these four species and those reported in Tribolium castaneum (63) could be divided into 12 clades and three major classes. The α-esterases significantly expanded in T. molitor, A. verucosus, and T. castaneum compared to dipteran and hymenopteran insects. In T. molitor, most COEs showed tissue and stage-specific but not a sex-biased expression. Our results provide insights into the diversity and evolutionary characteristics of COEs in tenebrionids, and lay a foundation for the functional characterization of COEs in the yellow mealworm.